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We study a new class of chaotic systems with dynamical localization, where gain/loss mechanisms 
break the Hermiticity, while allowing for parity-time VT symmetry. For a value "fvT of the gain/loss 
parameter the spectrum undergoes a spontaneous phase transition from real (exact phase) to com- 
plex values (broken phase). We develop a one parameter scaling theory for ^vT, a-nd show that 
chaos assists the exact PT-phase. Our results have applications to the design of optical elements 
with "PT-symmetry. 



Introduction.— Systems exhibiting parity-time {VT) 
symmetry have been the subject of rather intense re- 
search activity during the last few years. This interest 
was motivated by various areas of physics, ranging from 
quantum field theories and mathematical physics [IHl] to 
solid state physics |5ji6j and classical optics [7HO]. A sur- 
prising result that was pointed out in some of these inves- 
tigations was the possibility that VT symmetric hamil- 
tonians H can have real spectrum, despite the fact that 
they can, in general, be non-Hermitian [3]. The depar- 
ture from Hermiticity, is due to the presence of various 
gain/loss mechanisms which occur in a balanced manner, 
so that the net loss or gain of "particles" is zero. Further- 
more, as some gain/loss parameter 7 that controls the de- 
gree of non-hermiticity of H changes, a spontaneous VT 
symmetry breaking occurs. At this point, 7 = 7737-, the 
eigenf unctions of H cease to be eigenfunctions of the VT- 
operator, despite the fact that H and the T'T-operator 
commute [3]. This happens because the T'T-operator is 
anti-linear, and thus the eigenstates of H may or may 
not be eigenstates of VT- As a consequence, in the bro- 
ken VT -symmetry phase the spectrum becomes partially 
or completely complex. The other limiting case where 
both H and VT share the same set of eigenvectors, cor- 
responds to the so-called exact VT -symmetric phase in 
which the spectrum is real. 

A promising realization of VT symmetric systems ap- 
pears in the frame of optics, where a medium with alter- 
nating regions of gain and loss can be synthesized, such 
that the (complex) refraction index satisfies the condi- 
tion n*{-x) = n{x) pfTT]. This kind of synthetic VT 
materials exhibits unique characteristics such as "double 
refraction" and non-reciprocal diffraction patterns, which 
may allow their use as a new generation of unidirectional 
optical couplers or left-right sensors of propagating light 
[3] . Recently, the interest in VT systems bursted further 
due to their experimental realization \V2\ I13| . In this 
respect, one of the emerging questions is how one can 
enhance the parameter regime for which exact T^T-phase 
is present, while at the same time provide a general the- 
oretical formalism for the behavior of 7-pT: in terms of 
system parameters like imperfections, system size, com- 



plexity of the underlying classical (ray) dynamics etc. 

In this Letter, we investigate the behavior of the exact 
T^T-phase in a new setting of systems, namely a class of 
Hamiltonians whose classical (ray) dynamics is chaotic 
while its quantum/wave analogue can show dynamical 
localization [T31 [H] , a dual phenomenon to Anderson lo- 
calization appearing in disordered media [16] . As a re- 
sult of this duality, our study (although performed in the 
framework of wave chaos systems) is directly relevant to 
disordered quasi-one dimensional systems like disordered 
arrays of optical fibers [171 IE]- We have developed a 
one parameter scaling theory for 7^7- and show that it is 
the only relevant parameter that controls the variation 
of ^-pT with N , the system size of a sample. Specifically, 

— ^^/3(^.p.^); where ^-pT = N-fvT (1) 

olog A* 

where /3 is a universal function of 7^7- alone. 

Furthermore, we have investigated the distribution of 
the T'T-parameter V{j-pt) in the localized and delocal- 
ized/chaotic regimes and found that it reflects the prop- 
erties of the respective system. Specifically, we have 
found that Vij-pr) is log-normal in the former case, 
while in the latter it follows a Wigner distribution, re- 
flecting the level repulsion characterizing systems with 
chaotic/diffusive dynamics. Our results have direct ap- 
plications not only to coupled optical T^T-elements but 
also to cold atoms moving in a complex T'T-potential 

m- 

Model.— The prototype model of quantum chaos is the 
celebrated Kicked Rotor (KR) which exhibits the phe- 
nomenon of dynamical localization (DL) [T11[T5]. Specif- 
ically, it has been shown that the quantum suppression 
of classical diffusion taking place in momentum space is 
a result of wave interference phenomena similar in nature 
to the ones responsible for Anderson localization in ran- 
dom media. We study a variation of the Kicked Rotor 
(the VTKR) [HIinHll], defined by the time-dependent 
Hamiltonian |23j 
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where are a pair of canonical variables and q G 

[— TT, tt) [24]. The kicks have strength Kq and period T, 
which we set to unity without loss of generality. In or- 
der to avoid any integrable regions in the classical phase 
space (at 7 = 0) we take Kq > 5. We mark that in 
the framework of geometric optics Eq. ^ describes the 
propagation of light ray along a chain of optical elements 
equally placed along the axis of propagation t, in distance 
T from one-another [211 [55]. Furthermore we assume that 
the elements are purely refractive and ideally thin with 
a variation only in one transverse direction q. The phase 
space variable p = nodq/dt is proportional to the slope 
of the ray, while no is the free-space refractive index. 

The wave (quantum) dynamics of this system is de- 
scribed by the one-period evolution operator 

U = exp{-ip^/4n) exp{-ikV{q)) exp{-ip^/4n) (3) 

where k — Ko/h, p ^ Til ^ —iTid/dq and —N/2 < 
I < N/2. For 7 = 0, it was found that the eigen- 
functions ipi = (llip) oi lA, are exponentially localized 
in momentum space with a localization length ^ = 
limAT^oo l/Ef IV'; I" ~ fcV2 [Mlin]. They are solutions 
of the eigenvalue problem 
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A — exp(— ie) 



(4) 



where the eigenvalues A are unimodular at 7 = 0, and 
the phases e are referred to as quasi-energies. In the case 
of ?i = 2-kM/N with M,N integers, Eq. ^ defines a dy- 
namical system on a torus. The localization properties 
of the eigenstates are determined by the scaling parame- 
ter A = N /(^■. if A ^ 1 the eigenstates are exponentially 
localized while if A ^ 1 they are ergodically spread over 
the momentum space. 

VT — breaking scenario— In the exact VT phase (i.e. 
7 < 7737-) all eigenvalues are restricted to the unit cir- 
cle, resulting in a real quasi-energy spectrum (see insets 
of Figjl]). We find that the mechanism for transition 
to the broken VT phase is a level crossing between the 
pair of eigenvalues which are closest on the unit circle for 
7 = 0. However, rather than splitting into the complex 
plane symmetrically around the real line as in the case of 
Hamiltonian systems [^ , these pairs split in a logarithmi- 
cally symmetric manner away from the unit circle (one 
to the interior, the other to the exterior). Considering 
the first breaking pair of levels as an isolated two level 
system [57] we find that the first branching is described 
by |A±P cx exjp{±2y/^^j- — 7^). This square root singu- 
larity near the bifurcation point is quite universal of an 
exceptional point and applies both for A ^ 1 and A ^ 1. 
This is further confirmed numerically in Fig. [T] where we 
present the spontaneous PT-symmetry breaking scenario 
for two representative cases associated with localized and 
delocalized/chaotic parameter values. 

Scaling theory for 7^7 — We consider first the limit- 
ing case A ^ 1 where dynamical localization is domi- 
nant. To clarify the picture we start from the Hermitian 



limit 7 = 0. Imagine for the moment that the kick- 
ing strength k is zero. Then all states are 5— like func- 
tions localized at various momenta ~N/2 < I < N/2. 
There is an exact degeneracy of multiplicity two be- 
tween the states localized at ±/o, i-e. symmetrically 
around I — 0. For fc 7^ this degeneracy is lifted. The 
eigenstates whose centers of localization are a distance 
d = 2^0 3> ^ ~ fc^ apart, form a quasi-degenerate pair 
of symmetric /antisymmetric states [H]. Each has two 
peaks, near the momenta i^o, and decays exponentially 
V'(0 ^ (l/V?) exp(— |±/o — Z|/C) away from them (double 
hump states). Thus, the eigenstates in a V -symmetric 
KR are organized into pairs (doublets) ordered by quasi- 
energy difference, di < 62 < ■ ■ ■ . The splitting between 
quasi-degenerate levels is Si^ ^ (1/^) cxp(— 2Zo/^) while 
the energy separation between consecutive doublets is 
much larger, of the order of the mean level spacing of 
the system, A ^ 1/N. Specifically, the smallest energy 
splitting Amin = Si ^ (1/^) exp(— iV/^) corresponds to 
states that where originally located at the extreme points 
of the momentum "lattice", i.e. Iq = ztN/A. 

As 7 is switched on the eigenstates of each pair will 
initially preserve their T'T-symmetric structure [HI I28j . 
At 7 = 7737- ~ Aniin = Si, the two levels associated 
with Si will cross, breaking the T^T-symmetry (see Fig. 1 
inset). As 7 > j-pj- these modes cease to be eigenstates of 
the T'T-operator. Instead, the weight of each is gradually 
shifted towards one of the localization centers [51 [5S]. For 
larger 7 the next doublet (with splitting ^2) will come 
into play, creating a second pair of complex eigenvalues 
for 7 ~ (see Fig. 1 inset), etc. 

Let us now consider the opposite limit of A <C 1, where 
the eigenstates are ergodically spread all over the sys- 
tem. In this case, the picture of doublets with exponen- 
tially small energy splittings is not valid and 7^7- be- 
comes of the order of the minimal level spacing, Amim 
in the corresponding Hermitian problem. This state- 
ment follows from perturbation theory with respect to 
7. The unperturbed (i.e. 7 = 0) energy levels are real, 
and are separated by intervals of order 1/N, so that 
Aniin — l/N. Finite 7 leads to level shifts proportional 
to 7^ (the first order correction vanishes due to VT- 
symnietry) and for 7 — j-pj- ~ Amin the perturbation 
theory breaks down, signaling level crossing and the ap- 
pearance of the first pair of complex eigenvalues. Thus, 
the energy scale for the T'T-threshold in the N/^ <C 1 
limit (7757- ~ 1 /N) widely differs from that for N/S, > 1 

(7pr^(l/0e-^/«). 

Combining both cases, we conclude that in the two 
limits of weak and strong localization we have that 



Ivr = /(A) 



0(1) for 
Acxp(-A) for 



A < 1 
A > 1 



(5) 



Our numerical results for the VTKT model, are re- 
ported in Fig. [1] and are in excellent agreement with 
the above theoretical predictions. Moreover, they clearly 
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show that the scahng function f{x) is regular and inter- 
polates smoothly between the two limiting cases. This 
allow us to conclude that 



(6) 



which can be rewritten in the form of Eq. Q. This is 
the main result of the present Letter, as it allow us to 
postulate the existence of a /3-function for the 7^7- of 
generic chaotic (or quasi-lZ^ disordered) systems. In the 
remainder of the Letter, we will focus on the statistical 
properties of ^-pT ^ & function of the parameter A. 
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FIG. 1: (color online) Scaling behavior of 7^7- for the VTKR 
defined by Eq. (|3| for various M, and iCo > 5 with = 63 
(pink triangles), 127 (red circles), 225 (green squares), and 
511 (blue diamonds). All data shows nice collapse to the the- 
oretical prediction in Eq. ([bJi (dashed black line) . Inset: Para- 
metric evolution of eigenvalue magnitudes for Chaos (D.L.) in 
the upper (lower) figure with N = 127, and A = 0.05 (10). 
The first branching pair is responsible for the transition to 
broken VT phase, and in both cases follows the predicted 
functional form (see text) shown with blue circles. 

Distribution of j-p-j-.— The above discussion pertain 
only to the behavior of a 'typical' system. A full the- 
ory however must be formulated in statistical terms and 
deal with probability distributions Vij-pr)- To this end 
we exploit the equivalence between 7737- and Amin which 
is confirmed numerically in the inset of Fig. 2a for A 
values spanning the whole interval from the localized to 
extended regime. Thus, we instead analyze the distribu- 
tion P(Ainin)- For better statistics, an ensemble of V- 
symmetric KR systems has been created by randomizing 
the phases of the kinetic part of the evolution operator 
given by Eq. ^. In all cases, the numerical distribution 
7'(Ainin) involved more that 10^ data points for statisti- 
cal processing. 

We start our analysis from the localized regime A ^ 1. 
There are several sources of fluctuations in 61: fluctua- 
tions in the position and energy of the relevant localized 
states, as well as what can be termed "fluctuations in 
the wave functions" . By this we mean that a localized 



wave function exhibits strong, log-normal fluctuations in 
its "tails", i.e. sufficiently far from its localization cen- 
ter [29]. This latter source of fluctuations appears to be 
the dominant one and it immediately yields a log-normal 
distribution for Si (see Fig. [2]), since Si is proportional to 
the overlap integral between a pair of widely separated 
and strongly localized states [53] . We confirm numer- 
ically that for increased A we do in fact approach such a 
distribution (see Fig. [2]) . 

In contrast, in the delocalized regime A <C 1, the struc- 
ture of the eigenfunctions is random jMj [21] and they are 
ergodically spread over the momentum space. Thus there 
is no quasi-degeneracy for such states; instead, there is a 
level repulsion due to strong eigenstate overlap. This re- 
sults in a Wignerian distribution for the minimum energy 
difference, i.e. 7'(Amin) = (7r/2) A^in exp(-7rA^i„/4), 
which is dramatically different from the one found in the 
locafized regime (see Fig. [2^). Since P(Anii,i) ^ Vij-pr) 

we conclude that in the case of wave chaos V{"f-pi ^ 

0) — i.e. there always exists a 7-interval for which we 
will have an exact VT -phase. This observation can be 
used as a new criterion of wave chaos. Our numerical 
results for A = 0.01 are shown in Fig. [SJd, and are in 
excellent agreement with the above theoretical consider- 
ations. 
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FIG. 2: (color online) (a): Distributions of Amin for localized 
eigenfunctions displaying a convergence to log-normal behav- 
ior. Centers are shifted for ease of comparison. Inset: Linear 
relation between Amin and 7pt over roughly 12 orders of mag- 
nitude. Parameters and symbols are the same as in Fig. [T] 
(b): Distribution of minimum level spacings in the chaotic 
regime in which Wignerian behavior (blue line) is observed. 
Inset: The integrated distribution /(Amin) = /q"^""" P{x)dx 
in a log-log plot. The best fit (blue line) has power two. In 
both cases. A'' = 127. 



Conclusions - In conclusion, we have studied a new 
class of quantum chaotic systems with dynamical local- 
ization that also possess a T'T-symmetry. These sys- 
tems are described by a non-Hermitian hamiltonian due 
to the existence of well-balanced gain/loss mechanisms 
and show a spontaneous T'T-symmetry breaking, i.e. a 
transition from a real to a complex spectrum, for some 
value 7-p7- of the gain/loss parameter 7. We have de- 
veloped a one parameter scaling theory for the rescaled 
critical gain/loss parameter 7737- = j-p'j-N, and conclude 
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that there is a universal /3-function that depends only on 
"f-PT itself, which controls the variation of j-pT with the 
system size. Furthermore, we have analyzed the distri- 
bution V{'j-pf) in the localized/delocalized regimes, and 
show that it drastically differs in these two limits. In the 
former case it is log-normal while in the latter it follows 
a Wigner statistics reflecting the chaoticity of the under- 
lying classical dynamics. Our study opens the way to 
quantify the spontaneous breaking of the PT-symmetry 
in terms of universal /3-function. 

The results presented here are based on a simple con- 
nection between 7737- and the (minimal) level spacing 
Ainin (see inset of Fig. [2^) which is inspired by an iso- 
lated two level matrix model (see section on the VT- 
breaking scenario). Although this is the most generic 
type of scenario, describing a large number of physical 
realizations associated with classically chaotic (or quasi- 
ID disordered) systems, there are many other interest- 
ing cases which needs to be explored. For example, in 
the T'T-symmetric KR model presented above, we used 
a discontinuous (at the boundaries) imaginary potential 
luiV{q) which gives a lower bound on the behaviour of 
"f-pj- (worst case scenario). However, preliminary analy- 
sis [28] shows that if lmV{q) is a continuous and analytic 
function the T^T-symmetry breaking scenario can be dif- 
ferent. In such a case the matrix elements of lmV{q) be- 
tween V- doublets are exponentially small in the regime 
of dynamical localization. Thus, according to the lowest 
order perturbation theory, the energy levels of these two 
states remain nearly parallel (as a function of 7), and 
hence typically 7pt ^ \nin = ■ We also note that one 
can observe spontaneous anti-breaking of PT-symmetry, 
i.e. for some 7 ^ 7pt, a pair of non-unimodular eigenval- 
ues recombines again into a pair of uni-modular eigenval- 
ues, and sometimes, (but more rarely as N is increasing) 
one may even find situations for which all levels simoun- 
taneously become uni-modular (global recovery of the ex- 
act T'T-phase) [55] • Understanding of such anti-breaking 
mechanism could be of significant interest for optics ap- 
plications. Because of lack of space these results will be 
discuss elsewhere [28| . 

It will be interesting to extend this line of study to 
higher dimensions (possibly 3D disordered systems with 
a metal-to- insulator phase transition). We expect that 
our study will be of interest not only for the optics 
community but also for the atomic physics community 
where complex optical potentials have been recently con- 
structed [TO] . 
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